Measuring velocity of sound with nuclear resonant inelastic x-ray scattering 
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Nuclear resonant inelastic x-ray scattering is used to measure the projected partial phonon density 
of states of materials. A relationship is derived between the low-energy part of this frequency 
distribution function and the sound velocity of materials. Our derivation is valid for harmonic 
solids with Debye-like low-frequency dynamics. This method of sound velocity determination is 
applied to elemental, composite, and impurity samples which are representative of a wide variety 
of both crystalline and noncrystalline materials. Advantages and limitations of this method are 
elucidated. 
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Mechanical properties form an important part of our 
understanding of condensed matter. In many areas of sci- 
ence, measurements of sound velocity are used to study 
materials of both natural occurence and artificial fabrica- 
tions. For example, in the field of geophysics, the sound 
velocity is the most direct information we have about the 
Earth's interior. The standard approach to learn about 
the composition and structure of the Earth's interior en- 
tails measurements of sound velocities of candidate com- 
pounds. The results are then compared to seismological 
data to exclude or confirm a particular compound. In the 
following, we will describe the use of nuclear resonant in- 
elastic x-ray scattering (NRIXS) to measure the velocity 
of sound. 



The NRIXS method was introduced to probe the lat- 
tice dynamics of materials by employing low-energy nu- 
clear resonancesi^i In NRIXS experiments, only signals 
from nuclear resonance absorption are monitored, and 
for this reason the extracted quantity is specific to the 
resonant isotope. This technique provides the phonon 
excitation spectrum as seen by the probe nuclei, '^''^"'^ and 
in most cases one can extract the partial vibrational fre- 
quency distribution, a function often referred to as the 
partial phonon density of states (PDOS). The NRIXS 
method has been applied to various materials, e.g., thin 
films and multilayers )SiL£ nanoparticlcs,^-^*' crystals with 
impurities organic molecules jiSiiiMii^ proteins jiSiii 
samples under high pressures and samples of geophys- 
ical interests Most of these samples are compounds, 
and, while the obtained PDOS gives only part of the lat- 
tice dynamics, the low-energy portion of the PDOS pro- 
vides the Debye sound velocity of the whole sample. We 
will now show that, due to universal features of acoustic 
modes of harmonic solids, the low-energy portion of the 
PDOS is related to the Debye sound velocity in a simple 
way. 



The normalized phonon density of states is defined by 



(1) 



where the energy eigenstates of lattice vibrations Ei are 
labeled by quantum number I, and TV is the total num- 
ber of atoms in the solid. In the harmonic lattice ap- 
proximation, a PDOS, which is more relevant to NRIXS 
experiments, is given by^ 
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where v enumerates resonant nuclei, N is the total num- 
ber of resonant nuclei, k is a unit vector in the inci- 
dent photon direction, and e'\ are phonon polarization 
vectors. Equation |(2Jl shows that the vibrational polar- 
izations are projected onto the incident photon direction 
and in particular the vibrational modes with polarization 
perpendicular to the direction of the incident photon do 
not contribute. In the case of a single crystal, the mea- 
sured vibrational properties become dependent on the in- 
cident photon direction and were called "projected, " 
whereas in cases of polycrystalline or isotropic samples, 
the measured spectrum is an average over all directions. 
For a crystal in which resonant nuclei occupy only equiv- 
alent lattice sites, the quantity that can be extracted 
from NRIXS experiments is exactly described by Eq. |(5J . 
When the resonant nuclei occupy different sites, what 
can be extracted is an approximation of Eq. ^ . The ap- 
proximation is based on an average of phonon spectra for 
these different lattice sites. The closure conditions of the 
phonon polarization vectors guarantee the normalization 
of D(i?,k), i.e., its integration over all phonon energies 
is one. The orthonormality and closure conditions are 
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where a, /? denote the spatial components. These con- 
ditions hold for any harmonic solid, and the polarization 
vectors have to be specified for every atom in the solid, in 
contrast to the crystal case, where they can be reduced 
to a much smaller set for atoms in one unit cell only. 
The deviation of D{E,'k) from i^(-E) is contained in the 
behavior of the phonon polarization vectors e'^ and can 
be expressed in terms of a modulating function, 



(5) 



which we now determine for low-energy vibration modes. 

In a crystal, the acoustic modes form three branches 
described by phonon momentum. For disordered solids, 
we still expect hydrodynamic modes on length scales that 
are large compared to length scales characterized by inho- 
mogeneities in the material. In the appendix, we discuss 
the properties of such hydrodynamic modes. Our results 
show that linearly dispersing plane-wave modes exist on 
long length scales, and that these modes are described 
by momentum q and branch number s. Energies and 
atomic displacements associated with these modes take 
the following form 
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where c^s is the sound velocity, Pqs is a normalized po- 
larization vector, and are the atomic positions. The 
normalization factor ckqs can be determined as follows. 
For a normal mode / with energy Ei, the atomic dis- 
placements are given by^ 
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where e^" are phonon polarization vectors introduced in 
Eq. 1(21 and is the mass of atom /i. A comparison with 
Eq. ((BJ under consideration of the normalization condi- 
tion Eq. (PJ results in aqs = Ti/ yj2 m Eqs, where m is 
the average atomic mass. For the low-energy, hydrody- 
namic modes, we can therefore write 



-qs 
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In the small q, low-energy regime, we rewrite Eq. ^ 
by replacing the summation over phonon modes with an 
integration, J^i ~^ ^/(27r)'^ / q^dqdfiq, and we sub- 
stitute Eqs. © and © into Eq. @ to obtain 
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where m is the mass of the nuclear resonant isotope, and 
the Debye velocity vjj is defined as an average over all 
sound velocities 
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We see that for small energies, the modulating function 
becomes energy independent. If we define a projected 
sound velocity similarly by 
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we obtain the simple relationship 
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For an isotropic sample, the sound velocity does not 
have directional dependence, and we can further simplify 
Eq. H10|l and Hll|l to identify = ve,. In the case of 
a polycrystalline sample, averaging over all nuclear res- 
onant sites in Eq. (0) is equivalent to averaging Eq. @ 
over all directions of k. Thus, in both cases, we have 
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Finally, for an isotropic or a polycrystalline sample, 
the low-energy regime Eq. ||2Jl becomes 
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where n is the density of atoms. Equation (|14|) has ap- 
peared in a similar form in a NRIXS study of myoglobin 
and related biological compounds^ Here we have given 
a derivation of Eqs. ((T^ and In addition, Eq. l(T^ 

shows the dependence of the modulation function on the 
photon direction for anisotropic samples. A "mean sound 
velocity" was defined in the context of NRIXS in a pre- 
vious study.^^ It is identical to for isotopically pure 
samples, for which the PDOS becomes total DOS. Here 
we have shown that the sound velocity of a sample in 
general is obtained by applying a correction factor, the 
cube root of the modulation factor as given in Eq. 
to the "mean sound velocity" of Ref. |22|. 

The PDOS of a variety of samples has been measured 
by NRIXS. Here we show three examples representing 
crystals (bcc iron, hematite) and materials with point 
defects (^^^Sng.oiPdo.gg). The crystalhne samples are 
translationally invariant, and thus for them the basis of 
modes is in fact given by momentum eigenstates with po- 
larization vectors which are known to obey Eq. Q explic- 
itlyi^ While the point-defect case is expressly not trans- 
lationally invariant, by using the lattice Green's function 
technique, it is still possible to find closed form expres- 
sions for the polarization vectors, with Eqs. (|13|l and H14|) 
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FIG. 1; The measured PDOS divided by energy squared. 
The size of symbols indicates the statistical error bar derived 
from signal counts. These samples were measured with reso- 
lutions of 0.85 meV for "^Sno.oiPdo.99, 0.6 meV for ^''Fe203, 
and 1 meV for bcc iron. 



again being found to follow. We shall describe the details 
of this specific calculation elsewhere)^ while noting here 
that the modulation to the low-energy part of host lattice 
DOS is found to depend only on the mass ratio even in 
the event of a force constant change at the defect site.-^ 

The low-energy region of the PDOS divided by energy 
squared is displayed for the three samples in Fig. Q Bcc 
iron 95% enriched in ^^Fe is an example of the limiting 
case where iV ~ TV, and NRIXS provides the total rather 
than the partial DOS. Hematite enriched with ^^Fe repre- 
sents a situation, in which the resonant nuclei form only 
a part of the unit cell. The ^^^Sno.oiPdo.99 sample ap- 
proaches another limiting case in which the nuclear reso- 
nant isotopes occupy only a very small portion of lattice 
sites. In this limit, N ^ N, the average atomic mass m 
in the above equations is well approximated by the mass 
of a host lattice atom, and the sound velocity obtained is 
that of the pure host, rather than of the host /impurity 
system. These samples thus represent a very broad range 
of nuclear resonant isotope concentrations. 

We extracted numerical values for the velocities of 
sound by averaging 'D{E)/E'^ obtained from measured 
data in the region from zero to 5 meV according to 
Eq. (|14|l . The results are tabulated in Table U where 
we also compare our values with sound velocities from 



TABLE I: Comparison of velocities of sound, together with 
the factor (m/m) in Eqs. 1131 and 1141 . 
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''The lower limits from ref. 
"The upper limits from ref. I30l 

''Calculated with bulk and shear moduli from ref. Isi 



other sources as explained in the footnotes of the table. 
For iron and palladium, the results from NRIXS mea- 
surements are within the range of sound velocities that 
were obtained by other means. In the case of hematite, 
our results give clearly a lower value for vjj than is ob- 
tained from the measured elastic constants. It is difficult 
for us to judge the reliability of values for vd obtained by 
other methods that usually do not measure vd directly 
but rely on post-experimental data averaging. Values 
obtained with the NRIXS method seem to be lower or 
on the lower end of ranges given by other authors. This 
might be due to the fact that the phonon spectrum up to 
5meV (corresponding to a frequency of 1.2 THz) is used 
to derive the sound velocity. At such high frequencies, 
the phonon dispersion may already be nonlinear, which 
would typically lead to a reduction in the value of the 
obtained vjj. Improvements in energy resolution could 
make smaller phonon energies accessible, which would 
potentially provide a better measure of sound velocity. 

We have relied critically on the assumption of Debye 
behavior at small phonon energies. In particular, we re- 
quire Debye behavior to extrapolate to typical sound fre- 
quencies (IC^Hz) from the THz range, which is accessi- 
ble to NRIXS. However, in Fig.^ we see deviations from 
Debye behavior, which would correspond to a horizon- 
tal line. Besides nonlincarities in the phonon dispersion 
as mentioned above, the removal of the elastic contri- 
butions to the NRIXS spectra can be a source of sys- 
tematic uncertaintiesi^iSS421 Also the NRIXS technique 
relies on harmonic behavior of the sample to extract the 
PDOS from the measured data.'^^'^^ All these possibili- 
ties may contribute to the deviations seen in Fig. ^ In 
future studies aimed at low- frequency dynamics, we sug- 
gest measuring the resolution function simultaneously by 
nuclear forward scattering to improve the reliability of 
peak substraction. In any case, a highly accurate res- 
olution function should be available. Besides the peak- 
subtraction procedure and the energy resolution (usually 
defined as FWHM), the shape of the resolution function 
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is also of major importance. The access to small ex- 
citation energies is greatly improved if the tails of the 
resolution function can be minimized by design of the 
x-ray monochromator.^^ Furthermore, difficulties in the 
peak removal usually become more severe with increas- 
ing ratio of elastic to inelastic scattering intensities. In 
many cases, the relative strength of the elastic scatter- 
ing is reduced by saturation effects in the sample. ^ This 
highly desirable effect is less pronounced for samples with 
low concentrations of the resonant isotope or in cases of 
generally weak inelastic scattering, e.g., at very low tem- 
peratures. In addition to the special requirements for 
data collection near the elastic peak, an accurate mea- 
surement of the entire spectrum is equally important to 
achieve an accurate normalization. 

We have shown that nuclear resonant isotopes can be 
used to measure the sound velocity of a solid in the con- 
text of the harmonic approximation and Debye-like low- 
frequency dynamics. Our approach is valid even for very 
low concentrations of the nuclear resonant isotope, and 
therefore the probing nuclei provide information about 
the host lattice. We expect that the presented method 
will have significant impact in the scientific area of high- 
pressure research and, in particular, in the field of geo- 
physics, where sound velocity is of great interest. Also 
the strategic placement of resonant nuclei in artificial 
structures may provide insight into local atomic motion, 
which at low energies is thought to influence the elec- 
tronic noise in nanostructure devices. Thus NRIXS opens 
another venue to measure sound velocities of solids and 
can complement other techniques or even supercede es- 
tablished methods in those cases where they become too 
demanding or even impossible. 

This work and use of the Advanced Photon Source are 
supported by the U.S. Department of Energy, Basic En- 
ergy Sciences, Office of Science, under contract No. W- 
31-109-ENG-38, and by the State of Ilhnois under HECA. 



Appendix 

Our study of long-wavelength vibrational excitations 
starts with the continuum version of the equations of 
motion of a set of harmonically bound atoms, which reads 

w^m^u^ = - ^$^''u'^ , (15) 

V 

where are the force constant matrices and u'' is the 
displacement vector of atom /i for a vibrational mode 
of energy hio. We introduce mass and atomic density 
functions /9(x) and ?7(x) by 

p(x) - ^m^<53(x-rj 

77(x) = Y.^\^~v;) , (16) 



where is the position of atom [i. The continuum ver- 
sion that substitutes Eq. H15|) is then given by 

tjV(x)u(x) = - ^ ?/(x)$(x-x')?y(x')u(x')d3x' , (17) 

where we made the usual assumption that the force con- 
stant matrix depends on coordinate differences only. The 
original intentions to study long-wavelength excitations 
are served best by introducing Fourier transforms, e.g., 
u(q) = i^r[u(x)], and eventually expanding for small 
values of momentum. The transformed Eq. (|17|l reads 

/ p(q-q')ii(q')rfV = -(^ X 

y rXq - q')*(q')?)(q' - q")u(q") d^g'd^g" (18) 

The quantities ?7(q) and p(q) are closely related to the 
structure function S'(q) that is typically obtained from 
x-ray or neutron diffraction experiments. For crystals, 
these functions are described by a series of very sharp 
peaks at values given by the reciprocal lattice vectors. 
Disordered or amorphous materials do not produce these 
sharp peaks with the exception of the q = maximum 
which is not related to spatial order. One can assume 
that some of the salient features of Eq. (fT^ are captured 
by retaining the q = maximum only, i.e., we approxi- 
mate 

p(q) = (2.)3^<53(q) 

ry(q) = (2^f^8\ci^ , (19) 

where m and V are the average mass and volume per 
atom. Equation H18|) simplifies to 

mc^2Q(q) ^ _^|,(q)ii(q) ^ (20) 

and describes a class of solutions that should be com- 
mon for all harmonic solids. The previous equation per- 
mits us to label the modes with the value of q and a 
branch index s originating from the tensor character of 
the force-constant matrix, i.e., a; = cjqs. We note that 
$(0) = follows from the invariance of the vibrational 
energy with respect to a displacement identical for all 
atoms. Assuming that the force-constant matrix falls off 
sufficiently fast with distance an expansion of Eq. (|20|) in 
powers of g = |q| is permissible. The first-order term 
in this expansion will also vanish under the reasonable 
assumption of inversion symmetry of the force-constant 
matrix. For hydrodynamic modes, i.e., g — > 0, we may 
therefore use 

and long-wavelength vibrational excitations are de- 
scribed by a modified Eq. 1)20(1 

mu;^,u(q) =g2 0(^)„(q) ^ (22) 
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We obtain linearly dispersing modes with energies 
hujqs = Cqs witli sound velocities Cqs that depend on 
the direction of q. The matrix Q is related to the elastic 
tensor of the material (see for example ref. !29t) and is 
given by 

\ /q=0 M 

(23) 

N is the number of atoms in the sample. The eigen- 
values of Q provide us with the sound velocities, and 
the corresponding eigenvectors will not depend on q. 



The displacement field for a mode qs is then given by 



<.=«q.p<j.e*'i"-. , (24) 

where Pq^ is a normalized eigenvector of Q describing 
the polarization of the mode, and ckq^ is an appropriately 
chosen normalization factor. Equation H24(l shows that, 
even for disordered materials, hydrodynamic modes are 
equivalent to plane-wave excitations of the atomic dis- 
placements. 
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